Asymptotic Pade-Approximant Methods and QCD Current Correlation 

Functions 



F. Chishtiet V. Elias^ 
Department of Applied Mathematics, The University of Western Ontario, 
London, Ontario N6A 5B7, Canada. 

T.G. Steele* 

Department of Physics and Engineering Physics and 
Saskatchewan Accelerator Laboratory 
University of Saskatchewan 
Saskatoon, Saskatchewan S7N 5C6, Canada. 

February 1, 2008 



Abstract 

Asymptotic Pade-approximant methods are utilized to estimate the leading-order unknown (i.e., not-yet- 
calculated) contributions to the perturbative expansions of two-current QCD correlation functions obtained from 
scalar-channel fermion and gluon currents, as well as from vector-channel fermion currents. Such contributions 
to the imaginary part of each correlator are polynomials of logarithms whose coefficients (other than the constant 
term within the polynomial) may be extracted from prior-order contributions by use of the renormalization-group 
(RG) equation appropriate for each correlator. We find surprisingly good agreement between asymptotic Pade- 
approximant predictions and RG-determinations of such coefficients for each correlation function considered, 
although such agreement is seen to diminish with increasing rif. The RG-determined coefficients we obtain are 
then utilized in conjunction with asymptotic Pade-approximant methods to predict the RG-inaccessible constant 
terms of the leading-order unknown contributions for all three correlators. The vector channel predictions lead 
to estimates for the 0(a*) contribution to R(s) = [a(e + e~ — ► hadrons)/cr(e + e _ — » /i + fi~)] for three, four, and 
five flavours. 
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1. Introduction 



In phcnomcnological applications of perturbative QCD, one often encounters series of the form S — 1 + Rix + 
R2X 2 + R3X 3 + where x is the expansion parameter a s /n, and where only the first two or three coefficients Ri 
are known. Since calculation of the higher-loop diagrams for subsequent values of Ri becomes progressively more 
difficult, and since the expansion parameter x is often large (as in QCD sum rule applications), it is important to 
have a reliable method of estimating subsequent radiative corrections to the series S - that is, of estimating at 
least the first not-yet-calculated coefficient Ri. 

In the present paper, we employ the asymptotic Pade-approximant approach (already developed to determine 
successfully renormalization-group functions within QCD [1,2,3] and scalar field theories [1,4]) in order to estimate 
the leading uncalculated contributions to the imaginary (absorptive) parts of scalar- and vector-current correlation 
functions within QCD. As all but nonlogarithmic terms within these contributions can be extracted by renormaliza- 
tion group methods [3] , the asymptotic Pade-approximant estimates for coefficients of powers of logarithms can be 
tested against renormalization-group predictions for these coefficients. We find an astonishing degree of accuracy 
in such predictions, particularly for the coefficients of the highest and next-to-highest powers of logarithms within 
the leading uncalculated contributions to correlators. 

In Section 2, we consider the imaginary part of the scalar fermion-current correlation function, which has 
been calculated [5] to three subleading orders in a s /ir. The fourth order term may be expressed as a degree 
four polynomial in ln(s/so), where so is understood to be the continuum-threshold parameter of QCD sum-rule 
integrals in the scalar channel. We find good agreement between asymptotic Pade-approximant predictions and 
renormalization-group determinations for the four coefficients of ln k (s/so) [k = {1,2,3,4}] within the [a s /7r] 4 
contribution to the correlator, although this agreement diminishes as flavour-number n/ increases from 3 to 6. 

In Section 3, we consider the imaginary part of the scalar gluon-current correlation function, which has been 
calculated [6] to two subleading orders in a s /n. The leading uncalculated term is a degree-three polynomial in 
ln(s/so). Asymptotic Pade-approximant predictions for all three coefficients of ln k (s/so) [k — {1,2,3}] are found 
for n/ < 4 to be in excellent agreement (< 10% relative error) with values determined from the renormalization 
group, and within 22% relative error for rif = {5, 6}. 

In Section 4, the absorptive portion of the vector fermion-current correlator is analyzed. Once again, the leading 
unknown contribution to this correlation function, which has been calculated to three subleading orders in a s /TT 
[7] , is estimated by asymptotic Pade-approximant methods and compared to renormalization-group determination 
of the coefficients of ln k (s/ So) [k = {1, 2, 3}]. The agreement is found to be excellent for k — 3 and 2, but not for 
k — 1, where an approximate factor-of-two discrepancy is seen to occur. This discrepancy is explored further by 
expressing the asymptotic Pade-approximant estimate for the leading unknown contribution to the vector correlator 
as an asymptotic series in the variable L = —ln(s/so): R4 = d^L 3 + d^L? + d\L + ... . The two leading series 
coefficients d^^ of this large-L expansion replicate exactly the renormalization-group determination of the two 
leading polynomial coefficients ^3,2 within the i? 4 contribution to the vector correlator, i.e., the coefficients of 
— ln 3 (s/so) and ln 2 (s/ sq). The coefficient d\ of —ln(s/so) obtained from this large-L expansion, corresponding to 
an alternative Pade estimate of the k — 1 coefficient, is found to differ substantially from the correct value. However, 
these two asymptotic Pade-approximant approaches for estimating the k — 1 coefficient, as delineated above, are 
seen to straddle the correct result, suggesting for the vector correlator an overall insensitivity of asymptotic Pade- 
approximant methods to the k = 1 coefficient. A least-squares determination of the coefficients dk confirms this 
insensitivity to d\ . 

Finally, in Section 5 we utilize asymptotic Pade-approximant methods, in conjunction with renormalization- 
group determinations of ln k (s/s ) polynomial coefficients, in order to estimate those constant (i.e., k — 0) terms 
which cannot at present be extracted by renormalization-group methods. Such terms are of phenomenological 
interest for estimating higher-loop effects in QCD Laplace and finite-energy sum-rules. Moreover, such predictions 
can be tested against subsequent perturbative calculations. Of particular interest are predictions obtained for the 
O(af) term in R(s) = a(e + e~ — > hadrons) / a(e + e~ — > /i+/i~). 
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2. The Scalar Correlator for Fermion Currents 

We consider first the scalar current correlation function 

U s (p 2 ) = i J d 4 ye^y < 0\Tj.(y)j.(p)\0 >, (2.1) 
based upon the scalar fermion current 

j s (y) = VWO/) (2.2) 

The absorptive portion of this correlator may be expressed in terms of the ratio w = s/sq, where s = p 2 , and sq is the 
continuum threshold usual to QCD sum-rule analysis (i.e., the kinematic threshold above which purely-perturbative 
QCD alone is adequate to describe the correlation function): 



1 T TT 3S 

-ImU s = — 
7r 8ir z 



1 , a s (s ) . . / a s (s )\ 2 D , v 
H i?i(w)+ Ri\w) 

TT \ TT J 



+ (^) 3 i? 3 H + (^) 4 i? 4 H + ...j (2.3) 

The coefficients Ri(w), R 2 (w) and R 3 (w) have been determined for arbitrary flavour number [5]. For example, 
when rif = 3, these are given by 

Ri{w) = 17/3- 2 ln(w) (2.4) 

R 2 (w) = 31.8640 - (95/3)lnw + (17/4) (/n w) 2 (2.5) 

i2 3 (io) = 89.1564- 297.596 (229/2)(Zn w) 2 - (221/24)(/n w;) 3 . (2.6) 

This information is sufficient in itself to generate an asymptotic Pade-approximant prediction for i? 4 (u;) [3,4]: 



Rl(w) [R%(w) + R 1 (w)R 2 (w)R 3 (w) - 2Rf(w)R 3 (w)] 
R 2 {w) [2Rl{w) - R\{w)R 3 {w) - R 2 {w)R%(w)] 



This prediction, in turn, may be utilized to generate (via numerical integration) the first five (nonsingular) finite- 
energy-sum-rule moments [4] 



N k = (k + 2) [ dww K+l Ri{w) (2.8) 
Jo 

of the 0(a 4 ) contribution to the scalar correlation function. Substituting (2.7) into (2.8), we find that 

7V_i = 7544.9, iV = 2059.4, N x = 1158.4, N 2 = 833.47, iV 3 = 673.29 (2.9) 

The significance of these moments is that they may be used to estimate the polynomial coefficients of ln(w) in 
R±(w), which is known to be fourth-order in In w: 

i? 4 (u>) = do — d\ In w + d 2 (ln w) 2 — d 3 (ln w) 3 + di(ln w) 4 ; (2.10) 

Substituting (2.10) into (2.8) we see that 

7V_i = d a + di + 2d 2 + 6d 3 + 24d 4 (2.11) 

N Q =d Q + d x /2 + d 2 /2 + 3d 3 /4 + 3rf 4 /2 (2.12) 

iVi = d a + di/3 + 2d 2 /9 + 2d 3 /9 + 8d 4 /27 (2.13) 

N 2 = do + di/4 + d 2 /8 + 3d 3 /32 + 3d 4 /32 (2.14) 



N 3 = d + di/5 + 2d 2 /25 + 6d 3 /125 + 24d 4 /625. (2.15) 

The numerical values (2.9) already obtained for these moments lead to five linear equations in the five unknowns 
g?o — di. Their solution is 

d a = 252.5, di = 1339, d 2 = 1695, d 3 = 345.7, d 4 = 20.38 (2.16) 
The exact values for the coefficients d\ — d 4 may be extracted [4] from the renormalization group equation (RG) 







ImU(L(s ), x); (2.17) 



d d 

x = a s /n; L(s ) = ln(s /s) = -Inw; ( s — = — ) (2-18) 



ds dL 

As is evident from (2.3) and (2.4-6), the correlation function is in the following form 
3s 

ImH(L,x) — — [1 + (a + a\L)x 

OTT 

+ (6 + b x L + b 2 L 2 )x 2 + (c + CxL + c 2 L 2 + c 3 L 3 )x 3 

+ (do + d 1 L + d 2 L 2 +d 3 L :i + d i L 4 )x 4 + ...] (2.19) 

where dj, &, and Cj are known, and where the di are unknown. Given the known /3-and7-function coefficients 

0(x) = ~{t3 x 2 + fax 3 + /3 2 x 4 + ...), (2.20) 

lm {x) = -x(l + 7 ix + l2 x 2 + l3 x 3 ), (2.21) 

one can obtain via (2.17) the following set of equations for the d^. 

di = 3/3 c + 2/?i6 + I3 2 a a + 2 73 + 2 72 a + 2 7l 6 + 2c (2.22) 

2d 2 = 3/3 ci + 2/?i6i + /3 2 oi + 2 72 ai + 27161 + 2ci (2.23) 

3d 3 = 3/3 c 2 + 2/?!& 2 + 2 7l 6 2 + 2c 2 (2.24) 

4d 4 = 3/3 c 3 + 2c 3 (2.25) 

For nj = 3 the values of aj, 6j, are given in eqs. (2.4-6): 

a = 5.66667, ai = 2, 6 = 31.8640, 61 = 31.6667, 6 2 = 4.25 

c = 89.1564, ci = 297.596, c 2 = 114.5, c 3 = 9.20833 (2.26) 

Corresponding values for /3j and 7 i are [8] 

/3 = 2.25, /3i = 4, (3 2 = 10.0599, 
71 = 3.79166, 72 = 12.42018, 73 = 44.2628. (2.27) 

One then obtains from eqs. (2.22-25) the following RG determination of the O(aj) coefficients d\ — d 4 : 

di = 1563.0, d 2 = 1583.6, d 3 = 356.04, d 4 = 20.143. (2.28) 

The agreement between these numbers and the asymptotic Pade predictions (2.16) is astonishing close, ranging 
from a 14% relative error in the prediction of di to a 1.2% error in the prediction of d 4 . This close agreement 
supports the two-parameter asymptotic error formula [3] 

_ Rn+2 - R N+2 A 

5N+2 = ~^7 2 -N + l + B (2 ' 29) 

used to derive the asymptotic Pade-approximant prediction (2.7) in rcf [4]. 

The above agreement is not peculiar to nf = 3. Table 3 tabulates corresponding results for the rif = 4,5 
and 6 scalar-correlation functions. It is evident from Table 3 that the accuracy of asymptotic Pade-approximant 
predictions for d\ — d 4 is actually better for nf = 4 and 5 than it is for nf — 3. A diminution of accuracy becomes 
evident only when n/ = 6. 
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3. The Scalar Correlator for Gluon Currents 

We now consider the QCD gluonic scalar current correlation function = < >• This correlation 
function is of the same form as (2.1) but with j s (y) replaced with the RG-invariant gluonic current jc(y) 

h(y) - 3G(V) = ^^G%{y)G a ^{y) (3.1) 

The absorptive portion of this correlation function can be extracted to order-af from a previous three-loop calcu- 
lation [6] of the correlation function < (G 2 ) 2 >. Using our previous notation we find that 

Imn G (L,x) = ^[fo+^x + fox 2 +l3 3 x 3 ] 2 Im<(G 2 ) 2 > 

7T p 

= — [1 + [ao + aiL)x 
+ (bo + hL + b 2 L 2 )x 2 

+ (c + Cl L + c 2 L 2 + c 3 L 3 )x 3 + ...] (3.2) 

where x = a s /ir and L = ln(so/s), as before, and where the known coefficients ao,ai,6o,6i, and 62 are tabulated 
in Table 4. These coefficients are sufficient in themselves to determine the aggregate coefficients of x and x 2 in the 
correlator: 

Ri(w) = ao — ai Inw, (3-3) 

B.2{w) =60-61 Inw + b 2 [ln(w)} 2 , (3.4) 
2x 2 s 2 

ImIl G = ^- [1 + xR 1 (w) +x 2 R 2 (w) +x 3 R 3 {w) + ...] . (3.5) 

Asymptotic Pade-approximant methods may be utilized to predict that [4] 

o , ^ 1Rl{w) 

R ^ - Rl(w) + R 1 (w)R 2 (w) ' ( 3 - 6 ) 
This function may then be employed in the integrands of moment integrals 



P k = (fc + 1) / dww k R 3 (w) (3.7) 
Jo 

for k — {0, 1, 2, 3}. The numerical values of Pk can then be used, as before, to determine explicitly the coefficients 
c o, 1.2,3 that characterize R 3 (w) as a degree-3 polynomial in ln(w): 

R 3 (w) = Co — cilnw + c 2 (lnw) 2 — c 3 (ln w) 3 , (3.8) 

^0 = c + ci + 2c 2 + 6c 3 (3.9) 

A = c + Ci/2 + c 2 /2 + 3c 3 /4 (3.10) 

P2 = c + c 1 /3 + 2c 2 /9 + 2c 3 /9 (3.11) 

^3 = c + c 1 /4 + c 2 /8 + 3c 3 /32 (3.12) 

In Table 4, the values of these integrals Pk are tabulated for nj = {0, 2, 3, 4, 5, 6}, as extracted via (3.6) from (3.3) 
and (3.4). These integrals, in turn are sufficient to predict the coefficients {co, ci, c 2 , c 3 } by explicit solution of 
the set of linear equations (3.9 - 12). These predicted values of c _ 3 are labelled as asymptotic Pade-approximant 
predictions (APAP) and tabulated at the bottom of Table 4. 
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As in the previous section, all but one (co) of these coefficients may be extracted via the renormalization-group 
equation satisfied by the gluonic scalar-current correlation function 



— +(3(x)—)ImIl < ;[L..r) = 0. 



(3.13) 



where, as before, d/dL = so d/dso, and where (3(x) is given by (2.20) with the coefficients listed in Table 2. Upon 
application of (3.13) to the explicit form (3.2) for Imlla, one finds that 



ci = 2/3 2 + 3/?ia + Wo 
2c 2 = 3/3iai + 4/3 &i 
3c 3 = 4/3 fo 2 . 



(3.14) 
(3.15) 
(3.16) 



We have tabulated values of ci,c 2 ,C3 extracted via these equations as Ci(RGE), c 2 (RGE) and c^{RGE) towards 
the bottom of Table 4. In comparing asymptotic Pade-approximant predictions (APAP) for ci_ 3 to their actual 
values, as determined from RG (3.13), one finds accuracy within 4% for c\ and c 2 , and within 6% for c 3 when 
rif = 0. This accuracy diminishes somewhat as n/ increases, although the agreement remains striking even out to 
six flavours. Table 4's least accurate asymptotic Pade-approximant prediction, that of ci when n/ = 6, differs from 
the true (RG) value by only a 22% relative error. This agreement may be understood to confirm the applicability 
of the one-parameter simplification of the asymptotic error formula (2.29) 



Sn+2 — — 



N+l 



(3.17) 



obtained in ref. [1] and utilized in Section 5 of ref. [4] to obtain eq. (3.6). 

4. The Vector Correlator for Fermion Currents 

The vector-current correlation functions may be extracted from the Adler- function ( "D-function" ) presented in ref 
[7]: 



n„(Q 2 ) - 



(l + x) 



+ x< 



2\ \ 2 



2 ( /o 2 
+ B 2 [ln 



+ x° 
+ O^)} 
where x = a s /ir, as before, and where [7] 
A a = 1.98571 - 0.11529571/ 
A x = -1.375 + 0.083333371/ 

B a = -1.23954(^Q/) 2 /[3^Q^] + 18.2427 - 4.21585n/ + 0.0862069n^ 
/ / 

B x = -8.64820 + 1.04385n/ - 0.0192159nJ. 
B 2 = 2.52083 - 0.305556n/ + 0.00925926n^ 



(4.1) 

(4.2) 
(4.3) 
(4.4) 

(4.5) 
(4.6) 
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If we identify /j, 2 with the continuum threshold s , we find that the absorptive portion of this correlator is given by 
-ImU v = ~ V<2 2 f [1 + xR 1 (w) + x 2 R 2 (w) + x 3 R 3 (w)+x 4 R 4 (w)...] (4.7) 

f 

where w = s/so and where 

Ri{w) = 1 (4.8) 

R 2 {w) = A Q - {-2Ai)lnw (4.9) 

R 3 (w) = (B - n 2 B 2 ) - {-2B 1 )lnw + 3B 2 (lnw) 2 . (4.10) 

With this information, it is possible to repeat the calculation of Section 2 in order to estimate the 0(x A ) coefficients 
within Imll v . One can utilize (4.8), (4.9) and (4.10) within (2.7) to obtain an asymptotic Pade-approximant 
prediction of the function R 4 (w). Using this function, one can explicitly evaluate the moment integrals N-i, No, 
Ni, and N 2 , as defined by eq. (2.8). These four moments are sufficient to provide an estimate of the polynomial 
coefficients di within R 4 {w), which must be degree-3 in Inw: 

Ri(w) = d - dilnw + d 2 (lnw) 2 - d 3 (lnw) 3 . (4-11) 

The four linear equations relating moment integrals Nk to coefficients di are given by (2.11-14), with the constant 
di taken to be zero [the corresponding x A term for the scalar correlator was degree-4 in In w] . These four equations 
can then be solved to obtain asymptotic Pade-approximant (APAP) estimates of do, d\, d 2 , and d 3 . 

We have tabulated these estimates for n/ = {2,3,4,5} in Table 5. These estimated values for di- 3 can be 
tested against renormalization group determinations of these coefficients. The vector-current correlation function 
(4.7) can be parametrized similar to Ug(L,x) in (3.2): 



-ImU v (L,x) = -^V Q 2 f [l + x + (b + hLjx 2 + (c + CiL + c 2 L 2 )x 3 
3 / 

+ (do + diL + d^ 2 + d 3 L 3 )x 4 ...] (4.12) 

where L = ln(s /s)[— —Inw], and where the constants in (4.12) are related to those in (4.1) by 

bo = Ao, h = -2A U c = B - tt 2 B 2 , c x = -2B U c 2 = 3B 2 . (4.13) 

The correlator (4.12) has the same RG-invariance as (3.2), and is therefore a solution of (3.13). This equation 
predicts the following values for d\ , d 2 , and d 3 : 

di =3/3 c + 2/3ifeo + /32 (4.14) 

d 2 - (3#>ci + 2/3i&i)/2 (4.15) 

d 3 = (3 c 2 . (4.16) 

On can further show via (3.13) [or explicitly from substituting (4.2-5) into (4.13)] that 

h = fa, ci = 2/? &o + 0i, c 2 = 2 (4.17) 

in which case 

d 2 = 'S0 2 bo + ^0o0i, d 3 = 3 . (4.18) 

Table 5 displays RG determinations of di, d 2 , d 3 immediately below their APAP estimates. Striking agreement 
is evident from the Table between APAP and RG estimates of d 3 and d 2 for rif = {2,3}, an agreement which 
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deteriorates as nj increases. However, the APAP estimates of di seem to be a factor of two too large across the 
entire set of flavours considered. 

This discrepancy, however, is more indicative of insensitivity of APAP methods to the value of di and do 
for the vector correlator case, rather than of any overall failure of APAP methodology, as epitomized by eq. 
(2.7). Indeed, if we use eq. (2.7) directly (as opposed to the moment integrals Nk) to predict d 3 , d 2 , d\ and 
do, we not only obtain much different values for d\, but also results identical to the RG determinations of di 
and d 3 . To see this, consider the large-L asymptotic expansion of (2.7), with R1-4 as parametrized in (4.12) 
[-Ri = 1, i? 2 = b + biL, R 3 = c + ciL + c 2 L 2 }: 



i?4 — 



+ 



c 2 2 {b 2 + c 2 ) 



26? 
c|(l- 26 ) 
46? 

d APAP 'L 3 - 



46? 



L 3 + 

C\C 2 

61 



d APAP 'h 2 



_ [c|(l + 36o 
26? 

-0(L) 
d APAP 'L + 



+ 



26? 



Using (4.17) within the first two terms listed in (4.19), we find that 



d 



APAP' 



d, 



APAP' 



= 3/326 + 5/3 /3i/2. 



(4.19) 

(4.20) 
(4.21) 



These APAP' predictions, based entirely upon (2.7), coincide exactly with the RG determinations (4.18) of d 3 
and d 2 . In Table 6 we tabulate the corresponding predictions d APAP , as obtained via the large-L expansion (4.19), 
for the coefficient d\. Also tabulated are the prior (APAP) predictions of d x obtained via use of (2.7) within the 
integrands of moment integrals N k , along with renormalization-group (RG) determinations of d\. As is evident 
from the Table, the true (RG) value of d\ is straddled by the two Pade predictions, demonstrating the insensitivity 
of Pade methodology to d\, the "sub-subleading" 0(x 4 ) coefficient in the vector channel. 

This insensitivity of the di prediction is evident from a least-squares approach which finds values for dk which 
minimize the quantity 



1 

(d ,d 1 ,d 2 ,d 3 ) = J [R 4 {w) - (d - d\ \og(w) +d 2 \og 2 {w) - d 3 log 3 (w))] 2 dw 



(4.22) 



where R^w) is obtained by substituting (4.8-10) into (2.7). For n/ 
becomes the following quadratic form: 



3, the x 2 integral (numerically evaluated) 



X 2 (d ,di,d 2 ,d 3 ) 



227586.642 + 720^ + 12d 3 d + 2Ad 2 2 + 48eM 3 + 2d\ + 240d 2 d 3 + a% + I2d x d 2 + 2did 
+M 2 d - 237.6405798cio - 908.3757480(ii - 4402.296530(ii - 25479.03758^ (4.23) 



Minimization of the x 2 is equivalent to minimization of the matrix quadratic form 



(Ax - by 



where 



/ d \ 




( 


2 


2 


4 


12 \ 


di 






2 


4 


12 


48 


d 2 




4 


12 


48 


240 


V d 3 J 




\ 


12 


48 


240 


1440 J 



The singular value decomposition of A is 



A = 



S = Gib;. 



(To = 1481.970038 , a x = 10.52846941 , a 2 
( -0.008483667052 -0.2757146361 
-0.03336418987 -0.4565091835 
-0.1652682630 -0.8311862313 
\ -0.9856476345 0.1571949543 



U = V 



( 237.6405798 \ 
908.3757480 
4402.296530 
\ 25479.03758 / 



= 1.395692455 
0.8783851152 
0.2371489578 

-0.4115179052 
0.05341324608 



cr 3 = 0.1058005879 
-0.3903191722 ^ 
0.8568818869 
-0.3353637783 
0.03058614050 / 



(4.24) 



(4.25) 



(4.26) 
(4.27) 

(4.28) 
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and the quadratic form becomes 



3 



{Ax - bf = (Ey - b'f = £ (viVi ~ b' t f 



b' = 



U T b 



y = V T x 



(4.29) 



Although the quadratic form (and x 2 ) is minimized by 




2/o = -17.4587 



yi = -12.7421 



2/2 = -19.0195 



2/3 = -108.2255 



(4.30) 



the wide range of singular values cr ^> cr 2 > 03 implies that the x 2 depends strongly on 2/0 and 2/1 , but is relatively 
insensitive to 2/2 and 2/3- This broad x 2 minimum in the 2/2 and 2/3 directions compared with the sharp minimum 
in the 2/0 and 2/1 directions implies a comparatively large uncertainty in the extraction of 2/2 and 2/3 from the Pade 
approximation. Since 2/3 is both large and uncertain, the relation x = Vy implies that do and d\ (and to a lesser 
extent d 2 ) are dominated by the values of 2/2 and 2/3 rather than 2/0 and 2/1, and hence do, d\ are poorly determined. 
However, the 2/2, 2/3 dependence of the RG accessible dk can be eliminated to find the single linear combination 
independent of do that is well determined by the x 2 minimization: 



in extremely close agreement (0.33% relative error) with the Pade prediction (4.31). 

The insensitivity of Pade predictions of d\ for the vector correlator is to be contrasted with corresponding 
predictions for the scalar fermionic current correlation function considered in Section 2. If we apply the large-L 
asymptotic expansion of (2.7) to this scalar correlator, with i?i_3 as given in (2.4-6), we can predict APAP' values 
for cLt, c?3, d 2 and d\. These values are tabulated in Table 7 alongside the APAP estimates already listed in Table 
3. As is evident from the Table, both Pade methods predict quite similar values d&, d 3 and the sub-subleading 
coefficient d 2 , with the asymptotic large-L expansion (APAP') approach showing even greater accuracy than the 
moment-integral (APAP) approach delineated in Section 2. For the scalar correlator, the APAP' approach does 
not show signs of breaking down until the third subleading order of (2.7)'s large-L expansion, the coefficient d\, for 
which APAP' is considerably less reliable than APAP. This observation is confirmed by a x 2 minimization analysis 
of the scalar correlation function, which shows that the dk are much less sensitive to any poorly determined 2/i- 

Thus we see that different asymptotic Pade estimates do an excellent job of predicting leading and subleading 
0(x 4 ) coefficients for both the vector {d 2 ,d 3 ) and the scalar (d 3 ,d±) fermionic current correlators. However, the 
approaches diverge drastically in predicting the sub-sublcading term (dj) in the vector correlator, indicative of the 
limitations of the Pade method for this case. By contrast, the sub-subleading term {d 2 ) in the scalar coordinator is 
quite well predicted using either method, suggesting greater reliability of APAP predictions in the scalar channel 
for subsequent subleading terms. 



d 3 + 0.01950310416 d 1 + 0.1410349007 d 2 = -1.009606933t/ + 0.0310653405 2/1 = 17.23056 



(4.31) 



Using the RG values we find 



d 3 {RG) + 0.01950310416 d^RG) + 0.1410349007 d 2 {RG) = 17.17371 



(4.32) 
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5. Pade-Estimates of RG-Inaccessible Coefficients 



We have seen that the renormalization group equation may be utilized to determine all but one of the next-order 
coefficients in current correlation functions. Specifically, the coefficient do in scalar (2.19) and vector (4.12) fermion- 
current correlation functions, as well as the coefficient c in the scalar gluon-current correlation function (3.2), are 
not subject to RG constraints. However, knowledge of these constants is vital to a number of phcnomenological 
applications, such as higher order perturbative contributions to QCD sum-rules, which may be large at low So [3], 
or the O{o i s ) term in R = a(e + e~ — ► hadrons) / a{e + e~ — ► /i+/i~). 

Up until now, our emphasis has been on demonstrating how asymptotic Pade-approximant methods can estimate 
RG-acccssible coefficients within the correlators (2.17), of (3.2), and (4.12). However, it is evident from Tables 3,4,5 
and 6 that the accuracy of such estimates diminishes as the subscript index i of coefficients rfj and c, decreases. The 
APAP predictions for d in Table 5, for example, are quite suspect because of the insensitivity of APAP predictions 
to do and di in this channel, as discussed in the previous section. 

To obviate these difficulties, we can estimate do (cq) for fermion (gluon) correlators by first averaging the 
asymptotic Pade-approximant expression for R±(w) [R 3 {w)\ over the interval < w < 1, and by then removing 
the known contributions of RG-accessible coefficients <ii_4 {ci-3) to this average. For the scalar fermion-current 
correlation function, this procedure amounts to the use of RG-determined coefficients within (2.11): 



do= f R 4 (w)dw-d 1 (RGE)-2d 2 {RGE)-6d 3 (RGE)-2M 4 (RGE). (5.1) 
Jo 

The integral in (5.1) is just the integral iV_i tabulated in Table 3, where RG-valucs of di_4 are also listed. These 
new estimates [labelled do (APAP')] are presented in Table 8 for n/ values of phcnomenological interest, and are 
compared to the prior estimates of Table 3 obtained from higher-moment integrals. The two predictions are 
comparable for rif = 3,4, but quite far apart for rif — 6. 

Corresponding predictions for the scalar gluon-current correlation function may be obtained from rearrangement 
of (3.9): 

co= [ R 3 {w)dw-c 1 {RGE)-2c 2 {RGE)-6c 3 {RGE), (5.2) 
Jo 

with ci_ 3 (RG) and the integral (Pq) as given in Table 4. These new predictions for Co, labelled as co(APAP') in 
Table 9, are in good agreement with the previous Table 4 predictions co(APAP) for rif < 4, reflecting the somewhat 
better agreement of Table 4's predictions with RG values. 
Finally, do estimates for the vector fermionic correlator 

d Q = iV_a - d\ (RGE) - 2d 2 {RGE) - 6d 3 (RGE), (5.3) 

with all constants on the right hand side as tabulated in Table 5, are displayed in Table 10. Here the discrepancy 
with prior APAP estimates is quite large, but within a factor of 2 for the phenomenologically important rif = 5 
case. We reiterate, however that the Table 5 estimates of do are likely skewed by the factor-of-two disparity between 
corresponding estimates of d\ versus RG determinations of d\. Consequently, we regard the top line estimates of 
Tables 8-10 to be the ones to be tested against future calculation. 

As just one such example, the quantity R(s) is obtained entirely from the nonlogarithmic coefficients within 
(4.12):' 

R(s) = <r(e + e~ — > 7 — > hadrons)/a(e + e~ — ► /j + /i~) 

= Z^Q)[l + x + box 2 +c x 3 +d a x 4 ...}. (5.4) 
/ 

The coefficients 60 and c are given by (4.13) in terms of constants listed in (4.2-6), and are also tabulated in Table 
10. Thus, for five active flavours we predict that 

= ^I 1 + x + 1.40924a; 2 - 12.8046a; 3 + 31.5a; 4 ], (5.5) 

where the well-known first four terms [7,9] have been augmented by the rif = 5 APAP' estimate for do in Table 10. 
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This prediction for the x 4 term in (5.5) differs in both sign and magnitude from an earlier Pade-motivated 
prediction [10]. The prediction is also substantially less than the 73.5a; 4 term one would obtain applying (2.7) 
directly to the first three terms in (5.5). Such an approach, however, is based on the constants bo and Co only, 
corresponding to finding R±(w) only for the specific choice w = 1 within (2.7). By contrast, the result (5.5) devolves 
from an estimate incorporating our exact knowledge of all coefficients bi, c$, and di (except do) to find the average 
value of Ri(w) over the full range of w. 

In Table 10, the exact values for bo and Co are tabulated for nr — {2,3,4} as well. These values can be 
substituted into (5.4) to display rif = {3,4} expressions for R(s) which incorporate the APAP' estimate for do- 
lt is interesting to note that the predicted magnitude of the x 4 coefficients in R(s) is quite modest, suggesting 
that the accuracy of phenomenology based on the preceding three subleading orders of perturbation theory is not 
compromised by higher-order corrections. 
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n f =4 


n/ = 5 


n f = 6 


do 


17/3 


17/3 


17/3 


ai 


2 


2 


2 


bo 


30.5054 


29.14671 


27.7881 


h 


OA A A A A 

30.4444 


on noon 

29.2222 


28.0000 


b 2 


4.08333 


3.91667 


3.75000 


co 


65.1980 


41.7576 


18.8351 


Cl 


267.589 


238.381 


209.970 


C2 


104.384 


94.6759 


85.3750 


C3 


8.39352 


7.61574 


6.87500 



Table 1: Coefficients of 0(a s ), 0(a 2 s ) and 0(a%) in terms within the fermionic scalar-current correlation function, 
as defined in (2.1). 



nf = 


n/ = 2 


nf = 3 


nf = 4 


n/ = 5 


nf = 6 


11/4 


29/12 


9/4 


25/12 


23/12 


7/4 


51/8 


115/24 


4 


77/24 


29/12 


13/8 


2857 


48241 


3863 


21943 


9769 


65 


128 


3456 


384 


3456 


3456 


128 




3.93056 


3.79166 


3.65277 


3.51389 


3.37500 




14.8393 


12.4202 


9.94702 


7.41985 


4.83866 




61.8794 


44.2688 


27.3088 


11.0401 


-4.50240 



Table 2: Coefficients of QCD (i- and 7-functions, as defined in (2.20) and (2.21). 





n/ = 4 


nf = 5 


71/ = 6 




6217.2 


5050.4 


4045.6 


N 


1616.9 


1238.0 


920.15 


Nx 


874.58 


634.83 


437.31 


N 2 


609.96 


422.70 


270.86 


N 3 


480.56 


320.06 


191.81 


d (APAP) 


147.4 


64.18 


9.004 


di (APAP) 


1029.6 


744.9 


442.2 


di (RGE) 


1159.8 


791.52 


457.39 


d 2 (APAP) 


1398 


1177 


1067 


d 2 (RGE) 


1338.9 


1114.7 


910.31 


d 3 (APAP) 


307.4 


253.1 


174.3 


d 3 (RGE) 


305.7 


260.06 


218.82 


d 4 (APAP) 


16.63 


15.37 


17.29 


d 4 (RGE) 


17.312 


14.755 


12.461 



Table 3: Comparison of asymptotic Pade-approximant predictions (APAP) to renormalization-group equation 
(RGE) determinations of the 0{a, A s ) coefficients di — d 4 within the fermionic scalar-current correlation function. 
The integrals N i} as defined by (2.8) are obtained numerically from the APAP estimate (2.10) of R4, with the 
lower-order coefficients ai, bi, Ci (Table 1) determining i?i,2,3- APAP estimates of do — d 4 are obtained from these 
integrals via (2.11-15). RGE determinations of d 4 — d 4 are obtained via (2.22-25), with and 7, values as given 
in Table 2. 
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n f = rif = 2 n / = 3 n / = 4 n/ = 5 ft / = 6 



a 


1007 

44 


6919 

348 


659 
36 


4999 
300 


4123 
276 


367 
28 


ffli 


11 
2 


29 
6 


9 
2 


25 
6 


23 
6 


7 

2 


6o 


349.140 


246.434 


197.515 


150.210 


104.499 


60.3685 


bi 


3225 
16 


7379 
48 


2105 
16 


1769 
16 


4355 
48 


1153 
16 


&2 


363 
16 


841 
48 


243 
16 


625 
48 


529 
48 


147 
16 




11167.0 


6769.0 


4970.7 


3437.6 


2173.9 


1191.2 


Pi 


6966.7 


4052.1 


2870.2 


1875.5 


1075.7 


486.57 


^2 


5916.2 


3382.6 


2359.5 


1503.6 


823.73 


336.41 


^3 


5449.3 


3087.2 


2135.7 


1342.5 


716.80 


275.47 


c (APAP) 


4262 


2346 


1580 


950.4 


466.0 


144.5 


ci (APAP) 


4148 


2559 


1895 


1311 


805.8 


377.5 


ci (RG) 


4323.5 


2695.9 


2017.4 


1424.8 


915.12 


485.46 


c 2 (APAP) 


1145 


775.4 


622.8 


490.2 


376.3 


278.6 


c 2 (RG) 


1161.2 


777.76 


619.03 


480.73 


361.69 


260.75 


c 3 (APAP) 


77.92 


52.21 


41.71 


32.63 


24.92 


18.65 


c 3 (RG) 


83.19 


56.46 


45.56 


36.17 


28.16 


21.44 



Table 4: The first five rows list known terms of the gluonic scalar-current correlation function (3.2). The in- 
tegrals -Po-3j as defined by (3.7), are obtained numerically and utilized to obtain asymptotic Pade-aproximant 
(APAP) estimates of the next-order coefficients c -3 in the gluonic correlator. Renormalization-group equation 
(RG) determinations of c\, c 2 , and C3 are also tabulated to illustrate the accuracy of the APAP method. 
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n/ = 2 n/ = 3 n/ = 4 n/ = 5 



AT_1 


160.31 


118.82 


86.172 


63.940 




24.561 


20.699 


21.182 


26.106 


JVi 


11.817 


13.936 


19.911 


29.596 


iV 2 


10.173 


14.315 


22.373 


34.181 


do(APAP) 


20.51 


27.54 


40.22 


58.83 


di(APAP) 


-79.88 


-83.30 


-98.87 


-130.3 


d\ (RGE) 


-35.490 


-46.238 


-56.903 


-63.989 


d 2 (APAP) 


66.11 


51.93 


49.19 


62.01 


d 2 (RGE) 


59.701 


47.404 


36.561 


27.111 


d 3 (APAP) 


14.57 


11.78 


7.742 


1.899 


d 3 (RGE) 


14.114 


11.391 


9.0422 


7.041 



Table 5: Comparison of asymptotic Padc-approximant predictions (APAP) to renormalization-group equa- 
tion(RGE) determinations of the C(af ) coefficients di within the fermionic vector-current correlation function. 
The integrals Nk, as defined by (2.8), are obtained numerically from the APAP estimate (2.7) of R4. APAP 
estimates of do — c?3 are obtained from these integrals, as described in the text. 



71/ = 2 71/ = 3 71/ = 4 71/ = 5 

di(APAP') UA8 ^2420 ^33 -44.84 

di(RG) -35.490 -46.238 -56.903 -63.989 

di(APAP) -79.88 -83.30 -98.87 -130.3 

Table 6: Pade estimates cZi(APAP') for the vector correlator, as obtained via large-L asymptotic expansion of (2.7), 
are compared to the exact results di(RG) and to the Pade estimates di(APAP) of Table 5 obtained through use of 
(2.7) in integrands of moment integrals Nk (2.8). The two Pade estimates are seen to straddle the exact result. 
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n f = 3 rif = 4 n / = 5 n/ = 6 



d 4 (APAP') 


20.183 


17.323 


14.745 


12.434 


d 4 (RGE) 


20.143 


17.312 


14.756 


12.461 


d 4 (APAP) 


20.38 


16.63 


15.37 


17.29 


d 3 (APAP') 


355.5 


305.0 


259.2 


217.9 


d 3 (RGE) 


356.03 


305.73 


260.06 


218.82 


7 / A T~» A T~» \ 

d 3 (APAP) 


345.7 


307.4 


253.1 


174.3 




161 r 


1380 


1111 
1111 


916.2 


d 2 (RGE) 


1583.6 


1338.9 


1114.7 


910.31 


d 2 (APAP) 


1695 


1398 


1177 


1067 


di(APAP') 


1450 


244.7 


2582 


405.6 


d\ (RGE) 


1563.0 


1159.8 


791.52 


457.39 


d 4 (APAP) 


1339 


1030 


744.9 


442.2 



Table 7: Pade estimates d^APAP') of coefficients di for the fermionic scalar-current correlator, as obtained via 
large-L asymptotic expansion of (2.7), are compared to exact (RGE) results and to the APAP estimates of Table 
3. The APAP' estimates are seen to predict d 4 , d 3 , and d 2 with even better accuracy than the APAP estimates, 
but to suffer substantially diminished accuracy in predicting d\ . 



n f = 3 rif = 4 n/ = 5 n / = 6 
do (APAP') 195 130 115 156 
do(APAP) 252 147 64.2 9.00 

Table 8: Asymptotic Pade-approximant predictions of the do _coemc i cn t in the scalar fermion-current correlation 
function. do(APAP') is obtained from (5.1) using renormalization-group determinations of rfi_ 4 , as given in Table 
3. do (APAP) are previous predictions obtained in Table 3 without RG-inputs. 



rif = rif = 2 n/ = 3 n/ = 4 n/ = 5 n/ = 6 
c (APAP') 4021 2179 1442 8343 366! 55^6 
co(APAP) 4262 2345 1580 950.4 466.0 144.5 

Table 9: Asymptotic Pade-approximant predictions of the co-coefficient in the scalar gluon-current correlation 
function. co(APAP') is obtained from (5.2) using renormalization-group determinations of ci_3, as given in Table 
4. c (APAP) are previous predictions obtained in Table 4 without RG-inputs. 



rif = 2 rif = 3 rif = 4 rif = 5 

do (APAP') 1T2E) L90 15J 3L5 

c -9.14051 -10.2839 -11.6856 -12.8046 

b 1.75512 1.63982 1.52453 1.40924 

do(APAP) 20.5 27.5 40.2 58.8 

Table 10: Asymptotic Pade-approximant predictions of the do" coernc i en t in the vector fermion-current correlation 
function. d (APAP') is obtained from (5.3) using renormalization-group determinations of di_3, as given in Table 
5. do (APAP) are previous predictions obtained in Table 5 without RG-inputs. Exact values of bo and Co [7,9] are 
also displayed to facilitate use of (5.4) to obtain rif — 3,4 expressions for R(s). 
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